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Abstract 

The aim of this work is the derivation of two approximated expressions for the two dimensional 
Gaussian Q-function, Q(x,y;p). These expressions are highly accurate and are expressed in closed- 
form. Furthermore, their algebraic representation is relatively simple and therefore, convenient to handle 
both analytically and numerically. This feature is particularly useful for two reasons: firstly because it 
renders the derived expressions useful mathematical tools that can be utilized in numerous analytic 
performance evaluation studies in digital communications under fading; secondly because the two 
dimensional Gaussian Q-function is neither tabulated nor a built-in function in popular mathematical 
software packages such as Maple, Mathematica and Matlab. 


I. Introduction 


Special functions are invaluable mathematical tools in natural sciences and engineering. Their 
wide use in analytical performance evaluation studies in wireless communications often allows 
the derivation of closed-form relationships for vital performance measures such as channel 
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capacity and probability of error. In addition, the majority of special functions are built-in 
functions in popular mathematical software packages. As such, the algebraic representation of 
any related expressions as well as their computation have been significantly simplified. 

It is recalled that special functions play a significant role in the area of digital communi¬ 
cations [11]-[13] and [16]-[19] and the references therein. Such a special function is also the 
two dimensional Gaussian Q-function. It is denoted as Q(x,y;p ) and is constituted by the 
two-dimensional Gaussian integral, [l]-[5]. Alternative integral representations for Q{x, y: p) 
along with performance bounds were reported in [6]-[10]. Furthermore, an exact infinite series 
representation, which is expressed in terms of the incomplete gamma function and elementary 
functions, was derived in [11]-[12]. Nevertheless, it is noted that this series is the only explicit 
relationship for the two dimensional Gaussian Q-function since Q(x,y, p) is neither expressed 
in terms of other elementary and/or special functions, nor it is included as a built-in function 
in popular mathematical software packages such as Maple, Mathematica and Matlab. As a 
consequence, both its computation and analytical tractability are not as straightforward as in the 
case of most special functions in communication theory. 

Motivated by this, the aim of this work is the derivation of two tight approximations for the 
Q(x,y ; p) function. The high accuracy of the offered expressions is validated through compar¬ 
isons with results obtained from numerical integrations as well as by the identity Q(x , y\ p) = 
Q(x)Q(y). An important feature of the proposed approximations is that they are expressed in 
closed-form and have a rather simple algebraic representation. As a result, they are expected 
to be useful in various analytical studies relating to the performance evaluation of digital 
communications over fading environments. 


II. The Two dimensional Gaussian Q-function 


The two dimensional Gaussian Q-function is defined as: 



( 1 ) 


According to [2], the above relationship can be also represented in the more desired Craig form 
as follows, 
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Q(x,y, p) = 
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For the special case that x = y, equation (2) reduces to [2], 


Q(x,x-,p) 
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whereas for the case that p — 0, it reduces to the product of the one-dimensional Gaussian 
Q-functions of its arguments, namely, 


Q{x,y,p = 0) = Q(x)Q(y) 

An exact series representation for Q(x,y, p) was reported in [11]-[12], namely, 


(4) 


Q{x,y,p) 
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(5) 


III. Closed-form approximations for Q(x,y; p) 

It was mentioned above that the only explicit representation of Q(x, y: p) is the infinite series 
in (5). The algebraic form of this series is relatively simple. However, the fact that it has an 
infinite form typically raises convergence issues. More specifically, the number of terms required 
for truncating the series adequately is sufficiently related to its parameters. This is evident by 
the fact that for small values of x, y or p, only a few number of terms are required to truncate 
the series with small error. On the contrary, as as the value of y and p increases, the numbers 
of terms required for a small truncation error also increase. As a consequence, the required 
computation time increase as well. 

An effective mean to resolve such issues is the derivation accurate approximations that are 
expressed in closed-form and are also analytically tractable. Such approximated expression can 
be ultimately derived with the aid of two approximations for the one dimensional Q-function, 
Q{x), that were proposed in [11]-[13], namely, 


Q(x) — 0.49e' 


,ii _ 
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( 6 ) 
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and 


Q(x) ~ 0.208e- a87fe2 + OlSe" 0 - 525 * 2 + O.Me" 7 - 25 * 2 


(7) 


The behaviour and high accuracy of (6) and (7) with respect to Q(x) are illustrated in figure 



Fig. 1. Comparison of Q(x) in eqs. (6) and (7) with Q(x) in Matlab 


1. This can be also justified by figure 2 which depicts the level of the involved absolute error 
and absolute relative error, namely, 


e a =| Q(x) - Q{x ) 


( 8 ) 


and 
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ar 


Q(x) - Q{x) 
Q(x) 


(9) 
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Fig. 2. Absolute error, e a , and absolute relative error, e ar , between Q(x) and eqs. (6) and (7) 


A. A first approximation for Q{x. y: p) function 

A closed-form approximation for the Q(x, y; p) function can be derived as follows: The double 
integral in (1) can be alternatively expressed as, [11], 


Q(x,y,p) = 


2 Trv/r 


r ° 

r 00 u? puv 
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dv 


( 10 ) 


Consequently, by integrating once by part the inner integral and carrying out some straightforward 
algebraic manipulations, the following equivalent representation is deduced, 


Q(x,y, p) = 


e~*Q [ 1 dv 


( 11 ) 


\/27T J x \ y/l ~ P 2 

Thus, it is evident that an explicit expression for Q(x,y, p) is subject to analytic evaluation of 
the semi-infinite integral, 



d)dx 


where a,b,c,d E M. To this effect, by substituting (6) into (11), one obtains 


(12) 
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0.49 f°° v 2 8 y ~ pv „ (y-pv) 2 

Q(x,y ; p) ~ ._ / e 2 e 13 v^e 2 (i-^) dv (13) 

V 27T x 

which after algebraic manipulations yields the following closed-form expression, 



Fig. 3. Comparison of Q(x)Q(y) with Q(x,y\p = 0) in eq. 14 


where, 


0.49 8v v2 s 2 / 
Q(x,y,p)~—=e 13 Vi- p 2 e xvA- 

\M V 


A= 1 + 


5 

7 ^ 


i — p 2 


and 


(14) 


(15) 


5 = 


8p 


13>/r 


+ 


py 


p z 


p z 


(16) 


It is recalled here that the Q(x, y. p) function is not a built-in function in mathematical software 
packages such as Maple, Mathematica and Matlab. However, the following identity holds by 
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Fig. 4. Comparison of Q(x)Q(y) with Q(x,y,p = 0) in eq. 14 


definition, Q(x,y;p = 0) = Q(x)Q(y). As a consequence, the accuracy of equation (14) is 
assessed by means of the this identity in figures 3 and 4. One can clearly observe the excellent 
agreement between the plotted curves for almost all cases. This is also evident by the level of 
the involved absolute and absolute relative which for all cases is less than 5%. 


B. A second approximation for Q ( x, y: p ) function 

By following the same methodology, a second closed-form approximation for the two di¬ 
mensional Q-function can be derived. To this end, by making the necessary transformation of 
variables in (7) and then substituting in (11), the following relationship is deduced, 
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Evidently, the derivation of a closed-form expression for (15) is subject to analytic evaluation 
of the three involved integrals. Notably, these integrals have the same algebraic form as the 
(13). Therefore, after some straightforward algebraic manipulation, one obtains the following 
closed-form expression, 



Fig. 5. Comparison of Q{x)Q(y) with Q(x,y;p = 0) in eq. 16 


The behaviour of (16) is illustrated in figures 5 and 6 along with the Q(x ) function. One 
can observe the excellent agreements between the theoretical and approximated results. This 
agreement is also evident by the level of the involved absolute error and absolute relative error 
which is less than 4% for almost all parametric values. 
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where, 























Fig. 6. Comparison of Q(x)Q(y) with Q(x,y;p = 0) in eq. 16 
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IV. Closing Remarks 

In this work, we derived two novel closed-form approximations for the two dimensional 
Gaussian Q-function, Q(x, y. p). Two critical features of these expressions are their high accuracy 
and their relatively simple algebraic form. These features ultimately render them rather convenient 
to handle both analytically and numerically. The latter is particularly important since the two 
dimensional Gaussian Q-function is not a built-in function in popular mathematical software 
packages. As a consequence, the derived expression can be considered a useful mathematical 
tool in various analytical studies in natural sciences and engineering, in general, and particularly 
in applications in communication theory relating to digital communications over fading channels. 
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